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Abstract 

Gamow's approach to exponential decay of meta-stable particles via com- 
plex 'eigenvalues' (resonances) of a Hamiltonian is scrutinized. We explain 
the sense in which the non-square-integrable 'eigenfunctions' that belong to 
these resonances (Gamow functions) are relevant for the time-evolution of 
square-integrable wave functions. For concreteness we study a one dimen- 
sional square-well potential with a trapping region K and the evolution of 
wave functions, whose support is initially inside of K. It is shown that the 
sum over the first few time-evolved Gamow functions restricted to K yields 
an approximation for the evolution of these initial wave functions within the 
trapping region. The approximation is good for all times for which exponen- 
tial decay prevails. 
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CHAPTER 1 



Introduction 



The exponential law appearing in connection with the phenomenon of a-decay is experimentally 
well verified. Its quantum mechanical description has a long-standing history and remains the 
subject of research in mathematical physics to this day George Gamow [4] was the first to study 
a-decay within non-relativistic Quantum Mechanics in 1928. He considered solutions G to the 
stationary Schrodinger equation, which correspond to a complex 'eigenvalue' e = E — i\ rather 
than a real one. The Schrodinger evolution G e~ tet then immediately implies the exponential 
decay of the 'Gamow function' G with decay rate T. But due to the continuity equation this 
behavior is only possible when G increases exponentially with increasing x. Therefore, the Gamow 
function is not square-integrable and has no direct physical relevance. This raises the question 
whether there is a way to make sense of the Gamow functions. The following argument outlines 
our approach to the solution of this problem. 

The Schrodinger equation is a wave equation. Therefore, potentials whose core is separated from 
the outside by high barriers will confine wave packets for a long time. This is due to reflection 
of waves at barriers. Mathematically, this property of the potential becomes manifest in the 
existence of 'almost bound states' G or, as we have previously called them, Gamow functions. 
We already know that they solve 



by the lifetime V . 

On the other hand, the time evolution of a wave function ipQ within a potential V, which does 
not admit any bound states, can be determined by expanding tpQ in so-called generalized eigen- 
functions. This yields 



where the eigenfunctions <f>(k, x) are bounded solutions of the stationary Schrodinger equa- 
tion (1.1) for real e = k 2 /2 > 0. But if the 'eigenvalue' of the Gamow function is very close 
to the real axis, the generalized eigenfunction (j>(y2E, x) should mimic the shape of G(x), since 
both functions solve equation (1.1). This argument motivates that 



for all k in the vicinity of \J2E. 



Clearly, not every initial wave function i/jq will decay exponentially within a certain time regime. 
But Gamow's approach suggests that a square-integrable approximation of G yields the expo- 
nential law at least within a certain time interval. And a rough approximation is obtained by 







(1.3) 
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truncating all of G that lies beyond the core K of the potential, which yields 



4>o(k)= / x K G { x ) ( l ) {k,x)dx ~ c<q(k). 



Inserting this and relation (1.3) into the integral (1.2), the time evolution of the truncated Gamow 
function reduces to the Fourier transform of \7](k)\ 2 . Therefore, a considerable part of this work 
is devoted to the properties of the function r/. It will turn out that, in the vicinity of v2E, its 
squared modulus has Breit-Wigner shape 



Since the Fourier transform of the Breit-Wigner distribution is the exponential function, the time 
evolution of these initial wave functions is then approximately given by 



Thereby, we have not only found an initial wave function which decays exponentially within a 
certain time interval, but it also takes the shape of the Gamow function G. 

Now, this heuristic discussion needs to be turned into a rigorous one. In this regard, the first 
major step will be taken in Chapter 3, where the precise relation between the generalized eigen- 
functions <j) and the Gamow function G will be determined and where ij will be calculated. This 
will primarily rely upon the complex continuation of <p in the variable k. And the second major 
step will be taken in Chapter 4, where we will understand which initial wave functions tpo decay 
exponentially within a certain time regime. The basic technique utilized in the proofs of this 
chapter is the calculus of residues. 

Regarding the question whether Gamow functions have physical relevance, substantial progress 
was made by Skibsted with [9] and [10]. He also showed that the truncated Gamow function x^G 
decays exponentially within a certain time regime. And although expansions in eigenfunctions 
are used in his approach as well, his proofs do not rely on complex continuation. Nevertheless, 
he mentions this method as a possible alternative [9, p. 593], even though Skibsted [10, p. 47] 
considers complex continuation as difficult if rigorous results shall be obtained. However, the 
author believes that the expansion of e- iHt ip 

in eigenfunctions together with the application 
of the Residue Theorem to the resulting integral, is the most direct approach to a-decay. It 
clearly demonstrates how exponential decay arises and which role Gamow functions play in this 
regard. And in spite of the above concerns, we will obtain rigorous results using this method, at 
least for the specific potential studied in this work. In a sense, this thesis thereby explains the 
necessity of what was achieved by Skibsted [9]. Moreover, it will complement Skibsted 's work 
by explaining why truncated Gamow functions are a particularly good choice for proving the 
exponential decay of e~ tHt tpo within a certain time regime. 

Another question that arises in connection with Gamow's article is, in which sense a self-adjoint 
operator admits complex 'eigenvalues' or, as they are customarily called, resonances. In this 
respect, different definitions of resonances have been studied. The most popular ones are reviewed 
by Simon [8]. A recent review of Zworski [12] takes the widely accepted definition of resonances, 
as a pole of the resolvent (z— H)' 1 for granted. This allows for the study of asymptotic properties 
of, what could be called, the extended spectrum and the dependence of the resonances on the 
geometry of configuration space. 



\v(k)\ 2 



(k - Re V2e) 2 + (Im V2e) 2 





(1.4) 
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However, the main concern of this thesis is to understand the physical relevance of Gamow 
functions rather than questions regarding the extended spectrum. Therefore, a precise definition 
of Gamow functions will be provided in the next chapter. This will allow us to study their 
relation to physically relevant L 2 -functions in the framework of a one dimensional square-well 
potential. The motivation for choosing this potential was the article [5] of Garrido, Goldstein, 
Lukkarinen and Tumulka, whose main objective - contrary to this thesis - is to show that an 
initially localized wave packet leaves the square-well much slower than a classical particle. But 
before this specific potential is introduced and used to explain why an eigenfunction expansion 
serves to express e~ lHt ipo as an integral, we will show that square-integrable wave functions 
generically undergo exponential decay only within a certain time interval. In Chapter 3 we 
will then establish the precise relation between the generalized eigenfunctions and the Gamow 
functions for our specific potential. The function rj(k), which is calculated in this process, is 
of vital importance for Chapter 4. With its help it will be shown that not only truncated 
Gamow functions, but also bound states obtained by restricting the original Hamiltonian, decay 
exponentially within a certain time interval. The bound states will then allow us to extend these 
results to a general class of initial wave functions. And in this process we will understand why 
truncated Gamow functions are distinguished initial wave functions. The last chapter will then 
conclude with an outlook pointing to questions, which remain open. 



Notation and Units. The real part of a complex number z will be denoted by z' , its imag- 
inary part by z" and its complex conjugate by z. Although, the prime will denote the derivative 
with respect to a real variable from time to time, it will always be self-evident which meaning 
of the prime is referred to in the particular context. Moreover, 0(A") denotes a function /(A) 
which is asymptotically of the same order as A™ , that is 

r /(A) 
lim sup 

A^oo 

Apart from this, we choose units in which h - 



A" 



in 



< oo. 



1. 



CHAPTER 2 



Problem and Tools 



1. Statement of the Problem 

To approach the question whether the Gamow functions have physical relevance, we will now 
define them and rephrase the problem in mathematically precise terms. 

Gamow derived the exponential decay law from non-relativistic Quantum Theory, by considering 
functions ip that solve the stationary Schrodinger equation 

1 d 2 

-~^{x) + V{x)ilj{x) = e^{x) (2.1) 

for complex rather than real e. Since such a ip yields the solution 

V<i,-) = e- le V(-) (2.2) 

to the time dependent Schrodinger equation, it decays exponentially in time when Im e is negative. 
But due to the fact that (2.2) needs to respect the continuity equation 

d t \t/j\ 2 + div (Im tpVip) = 

as well, the flux associated to this solution must be outgoing such that both terms on the left 
hand side cancel. This motivates the following definition of Gamow functions. 

Definition 1. Given the one- dimensional stationary Schrodinger equation (2.1) with a poten- 
tial V having compact support and non-zero complex e having negative imaginary part. A func- 
tion G G C X (M) is called Gamow function, if it solves (2.1) and if there is a non-zero constant C 
such that 

G{x) = Ce +iz W Vxe(suppV) c , 

where e = \z 2 = E — i^. 

Thus, Gamow functions are defined by regarding the Schrodinger equation as ordinary differential 
equation on C 1 rather than as an operator identity on a Hilbert space. Since the notion of self- 
adjointness thereby loses its meaning, it should not be disturbing that the Gamow 'eigenvalue' e 
is complex. 

From Definition 1 it becomes evident that Gamow functions do not lie in the Hilbert space L 2 (M). 
For, the negative imaginary part of z directly implies that G{x) must increase exponentially as 
\x\ tends to oo. Therefore, the question whether the Gamow functions have physical relevance 
can be rephrased to whether there is a relation between Gamow functions and Hilbert space 
solutions of the time dependent Schrodinger equation. 

As already outlined in the introduction, we will find a relation between G and the generalized 
eigenfunctions 4>. And although these eigenfunctions do not belong to L 2 (M) either, they can be 
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used to express the time evolution of a square-integrable wave function ipo as the integral 

e~ im ip = J Mk)<f>(k,-)dk. 

And this will allow us to establish a relation between physical relevant Hilbert space solutions 
of the time dependent Schrodinger equation and Gamow functions. 



2. Properties of the Schrodinger Evolution 

In the current section we will see that e~ lHt ip can not decay exponentially for t > 1 as well 
as t <C 1. This will follow from two well-known properties of the survival probability, which is 
defined as 

P M ,{t) = \{^e- iHt ^)\ 2 . 

The proof presented here, is based upon the arguments given by Simon in [8]. And it relies 
heavily on the spectral theorem for unbounded self-adjoint operators (see e.g. [1, Chapter 2]), 
such as the Schrodinger operator. Therefore, we will briefly review the contents of this theorem. 

Casually speaking, the spectral theorem generalizes the fact that every hermitian matrix can be 
diagonalized, to operators on Hilbert spaces. To make this precise, let H denote the Schrodinger 
operator, H denote the Hilbert space on which it is defined and let cr(H) denote its spectrum. 
Then the spectral theorem guarantees the existence of a finite countably additive measure \jl 
on cr{H) x N and the existence of a unitary map 

J 7 :H —> L 2 (a(H) xN,/i). 

This map is such that J-HJ 7 ^ 1 = h, where h is the multiplication operator associated to the 
function h(n, s) = s with s being an element of o~(H) . Since h is the analogue of a diagonal matrix, 
this justifies the casual statement of the spectral theorem given above. However, T develops its 
full power only when it is applied to the operator f{H) with / being an arbitrary bounded 
measurable function. In this case we find analogously to the identity from above J- f(H)!F~ l = 
f{h), which will be used extensively in the following two proofs. The fact that J- maps onto a 
direct sum of several L 2 (o~(H), /i)-spaces is due to the devision of TL into cyclic subspaces, which is 
necessary to prove the spectral theorem. But on every cyclic subspace, H is unitarily equivalent 
to a multiplication with h(s). Therefore, we will neglect this detail in the sequel without causing 
any harm. 

Having reviewed the spectral theorem, we can now turn to the first of the above-mentioned 
properties of the survival probability P^,. It concerns the short time behavior of e~ lHt ip. 

Lemma 1. Let H be a self- adjoint operator and let ip be an element of its domain. Then P^{t) 
is differentiate everywhere and -P^(O) = 0. 

Proof. The differentiability of P^(t) follows from the differentiability of 

A(t) = ^,e~ iHt i;). 

Since H is self-adjoint, the spectral theorem is applicable. This directly implies 

A(t)= [ |^V(e)| 2 e--V(d6). (2.3) 



2.2 Properties of the Schrodinger Evolution 
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Using dominated convergence and the inequality of Cauchy-Schwarz, we thereby see that A is 
differentiable if 

/ e\^(e)\ 2 e-^ n{de) = \^,e~ im H^)\ < ||Y>|| \\H^\\ < oo. 

Since this is the case for all ip in the domain of H, the survival probability P^ (t) is differentiable 
for all t in R. 

Using the spectral theorem again, can be rewritten as 

|{^e-^}| 2 = {i^}{^e^), 

which implies that P^{t) = P^(—t) for all t in M and ip in the domain of H. Moreover, an 
application of the inequality of Cauchy-Schwarz proves that P^(t) < P^(0). From these facts 
and the differentiability of P^, we can conclude that P^(0) = 0. □ 

The differentiability of P^ and the vanishing derivative at zero, immediately imply that 

P^t) ± Ce~ vt 

for small t. So even if e~ lHt ip contains a contribution of the form e~ rt ip, there must be an 
additive remainder that accounts for the non-exponential behavior of P^ for small t. Apart from 
this restriction on the short time behavior of e~ lHt ip, there is another property of the survival 
probability having implications for the long time behavior. 

Lemma 2. Let H be a self-adjoint operator, whose spectrum is bounded from below. Then the 
only element ip in the domain of H that satisfies 

P^t) < Ce~ r W 

for some C,T > is the zero element. 



PROOF. Since H is self-adjoint and ip is an element of its domain, equation (2.3) can be 
utilized again such that 

A(t) = [ \^(e)\ 2 e^ fi(de) = [ g^e^ fi(de), 

Ja(H) JK 

where g is equal to \ J r ^\ 2 on the spectrum a(H) and zero otherwise. Therefore, A is the Fourier 
transform of a L 1 -function g, which is supported only on the spectrum of H. By inverting the 
Fourier transform, we find 

g(e)= I A{t)e^ii{dt). 

And if we assume that 

\A(t)\ < Ce-£l*l 

for some non-zero positive constants C and T, this formula together with the theorem of dom- 
inated convergence implies that g is analytic on the strip {eo + irj\ e o £ £ (— r/2,T/2)}. 
But the spectrum of H is bounded from below. So g(e) vanishes for all real e smaller than the 
infimum of <r(H). Therefore, g must vanish on the whole strip, which is only possible if ip equals 
zero. □ 
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That there is no exponential bound on P^, can only be due to the asymptotics of the survival 
probability and thereby due to the asymptotics of e~ tHt tp. It reflects the fact that the wave 
function ultimately decays polynomially instead of exponentially. This is not surprising, since it 
is in accord with the following heuristic scattering theoretical argument (for a detailed discussion 
see [2, Remark 15.2.4]). The time evolved wave function at the fixed position x is given by 

e -^V(x) = r 1/2 J (-j=) <j> dK, 

where n/y/i was substituted for k. In the limit of t tending to oo the right hand side multiplied 
by t 1 / 2 turns out to be constant, which can be seen heuristically by interchanging limit and 
integration. For f> 1 we therefore expect that 



The previous lemmas also show that the heuristic discussion in the introduction, which lead us 
to 

e- iH ty(aO ~ G(x)e-?l*l, 
can only be valid on intermediate time scales. However, for very small as well as very large times 
the error terms, hidden behind the ~ symbol, will not be negligible anymore. In Chapter 4 this 
heuristic statement will be turned into a precise one. 

Let us conclude this section with the following aside. If exponential decay can prevail on inter- 
mediate time scales, we could just choose to big enough and consider 

P i ,(t + t) = \(i>,e- iHt e- iHt °^\ 2 

in order to ensure exponential decay even for times t <C 1. This suggests that by choosing the 
initial time to advantageous the result of Lemma 1 can be somehow circumvented. But in fact 
resetting the initial time results in the quantity 

\(e' iHto ^,e- iHt e- iHto ip}\ 2 

as opposed to P^(to + t), since we really want to compare the time evolved wave function e~ lHt t^t 
to the initial wave function ipt . And this quantity is identical to P^(t), which is why it can not 
decay exponentially for times t <C 1 either. 



3. Setting and Tools 



The specific potential that provides the framework for this thesis is given by 

V X (X) = XX[-a,a]( X ) with «,A>0 (2.4) 
and is illustrated in Figure 2.1. Due to the reflection of waves at barriers, it confines initially 
localized wave packets for a long time if A S> 1. As was already mentioned in the introduc- 
tion, V\ will therefore admit resonances. Since it is explicitly solvable too, V\ is a good choice 
for approaching the question whether Gamow functions have physical relevance. 

Clearly, the function V\(x) defines a multiplication operator V\ acting on L 2 (M). It is not dif- 
ficult to see that this operator is symmetric and bounded. And in particular, V\ has relative 
bound zero with respect to the free Schrodinger operator Hq acting on L 2 (M). The Kato-Rellich 
theorem (see [1, Theorem 1.4.2]) therefore immediately implies that H\ = Hq + V\ is self-adjoint 
with domain H 2 (R) = V(Hq). And due to the self-adjointness, the spectral theorem guarantees 



2.4 Resonances and Gamow Functions 
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Figure 2.1. Plot of the potential V\. 



the existence of a map J- yielding the spectral representation of H\. In the previous section we 
already used this map, in order to diagonalize e~ lHt . Regarding the time evolution of the wave 
function ipo, this only yields 

In order to get any further we need to determine J- explicitly. Fortunately, for our simple 
potential this can be done. One possibility to achieve J 7 , is to determine the integral kernel of 
the resolvent (z — Hx) which rather directly yields the spectral measure dP e . This strategy is 
pursued in [11, Chapter 23] to which we refer for details. However, after doing the calculations 
one will find that the map T yielding the spectral representation of H\ can be defined as an 
integral operator 

Fip(k) = / il){x)<t>{k,x)dx (2.5) 
Jr 

with respect to the Lebesgue measure. And the integral kernel turns out to be an element of L°°, 
which solves 

1 d 2 k 2 

- -—,(j)(k,x) + V4>{k,x) = e(f>{k,x) for e= — >0 (2.6) 

2 dx z 2 

and is therefore called generalized eigenfunction. Although, bona fide eigenfunctions usually 
enter into the explicit expression of T as well, this is not the case for V\. For this particular 
potential the pure point spectrum is empty. Thus, there are no bona fide eigenfunctions and the 
spectrum is given by 

a{H x ) = [0,oo). 

Before we continue, notice the similarity of equation (2.5) with the Fourier transform. It has 
earned T the name generalized Fourier transform. Clearly, the generalized eigenfunctions <p 
define the spectral representation of the perturbed Hamiltonian H\ in complete analogy to the 
way the plane waves e tkx define the spectral representation of Hq. And it is often helpful to 
remember this analogy. 

4. Resonances and Gamow Functions 

In the present section we will determine the Gamow functions that belong to the potential V\. 
Thereby, we will not only have an explicit example for their generic definition. But in the 
process of calculating the Gamow functions, we will also determine the approximate location of 
the corresponding 'eigenvalues'. Although these eigenvalues will be called resonances from now 
on, it will only be seen in the next chapter that their defining equation coincides with the formula 
which defines the poles of the S-Matrix thereby justifying this name. 
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An ansatz for the solution of equation (2.1) with potential V\, which satisfies the boundary- 
conditions of Definition 1 is given by 



G(x) 



with 



Ae -izx 
Be izx 
De izx 



+ Ce 



, x < —a 
, \x\ < a 
,x>a 



Vz 2 - 2A. 



Whereas the complex square root is defined such that Re ^/w > for all w in C \ (— oo, 0]. So 
whenever Im w ^ we know that Im \Jw ^ 0. For G to be an element of C 1 as required from 
the definition of Gamow functions, its first derivative and itself have to be continuous at x equal 
to ±a. This will be the case when the coefficients satisfy the following homogeneous system of 
linear equations 



V 



e 

-ize 1 





-ize 
—e 1 ' 



-ize 



ize 
—e~ 



ize 






e izc 

ize 1 " 









B 




C 


) 





(2.7) 



Due to the fact that this system has non-trivial solutions only when the determinant of the 
matrix vanishes, the equation 

2 



AAaz 



z + z 



z — z 

needs to have solutions for Gamow functions to exist. 



(2.8) 



Lemma 3. Given the potential V\ with fixed width a > 0, there is a Gamow function corresponding 
to every integer n which satisfies 

2 



1 < n < -\ / 2a 2 \ - 1. 



7T 



(2.9) 



Moreover, if such an n is fixed the corresponding resonance is given by 



2A + 



n7r 
2a~ 



n 2 7r 2 
nir 



1/2 



n 2 7r 2 



2aV2A 



4a 3 V2A 
+ 0(A- 1 ) 



2A + 



n 2 7r 2 
~4aJ~ 



-1/2 



+ 0(A^ 3/2 ) and 



respectively. 



PROOF. Due to the fact that exp(— in2ir) = 1 for all n in Z, the identity (2.8) represents a 
whole family of equations numbered by the integer n. Using the definition of z, this family can 
be expressed as 

4 



e UaV2Xk-in2n = ^ + ^ ' R 2 + ^ j 



where k denotes z/y2\. And with the help of some simple manipulations we can even bring it 
into the fixed-point form 

7T 1 / 

F n (k) := n -j= — i — -j= ln(/t + y k 2 + 1) = k, 

2ay2\ av2A 



2.4 Resonances and Gamow Functions 
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where In denotes the principal branch of the complex logarithm. Therefore, the Banach Fixed 
Point Theorem (see [3, Chapter 9.2.1]) can be applied to each of these equations in order to find 
their solutions by iteration. 

This theorem is applicable whenever F n is contractive on a closed subset T of C and T is invariant 
under F n . Due to the fact that the derivative of F n satisfies 

■wi 1 1 1 1 

\F = < 



aV2\\l + K 2 \ 1/2 ~ aV2\\^-\^\ 2 \ 1/2 ' 
the map is contractive on T = {k G C||ft| < K} with K being a constant which satisfies 



0<K<J1 ' 



2a 2 X' 



The requirement of having an invariant subset will restrict K further. In order for |.F„(re)| < K 
to hold, we need 

\F n {k)\ < \F n (K)-F n (0)\ + \F n (0)\ < 1 p. +\ n \—^—<K. 

av 2a -L \k\ 1 ztiv 2\ 



But since every term is positive, this can only be the case for non-zero k, if 

II 71 " V 

n = < K. 

2aV2\ 



Together with the upper bound on K from above, we thereby conclude that 

1+1 — I <2a 2 X (2.10) 



2 

needs to hold true for such a K to exist. Thus, for fixed parameters a and A the Banach Fixed 
Point Theorem can only be applied to those fixed-point equations whose index satisfies (2.10). 
If the parameters are such that this inequality is fulfilled, n = will always be allowed. But 
the corresponding fixed-point equation has k = as its unique solution, which leads us to the 
trivial solution of the Schrodinger equation. Since this was excluded in the definition of Gamow 
functions, we require |n| > 1. 

To prove the second part of the Lemma let a > be fixed, assume that A is big enough for (2.9) to 
allow resonances and fix an n which satisfies this condition. Then we can iterate the corresponding 
fixed-point equation in order to find the location of the nth resonance. If this is done up to second 
order starting with k^ = 0, we end up with a rather cumbersome expression. Expanding this 
expression in powers of A -1 / 2 and keeping only the first non- vanishing order of real and imaginary 
part yields 

K n = n /HT " Z 7~2T + A 7 alld Z n = ~ K 2 ATT + A 

2aV2A 4a z A 2a 2a 2 v2A 



respectively. If we use the definition of z and expand again we also get 

-,2^2 X 1 / 2 „2^2 / ^,2^2 \ -V 2 



where we implicitly used that Im y/z < whenever Im z < 0. And in order for Im z n to be 
negative, n is restricted to positive integers. □ 
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To determine the Gamow functions that correspond to the resonances just found in the previ- 
ous lemma, it remains to solve the homogeneous system of linear equations (2.7). A series of 
straightforward calculations yields 



' 2z — iza „— iz(x+a) 



G(x) = B < 



, x < —a 



e izx + ^ e iza e -iz(x-a) ) | x | < a 



2z iza iz(x-a) 
z+z^ e 



, x > a. 



Evaluating this expression at z = z n immediately implies that 



z„ c iz n a c —iz n (x+a) x ^ _ a 

Zn~\~Z n 



G n (x) = < 



cos z n x 



z n c iz n a c izn(x-a) 

< Zn~\~Z n 



, \x\ < a 



, x > a 



for odd n, while 



G n {x) = < 



i sin z n x 



, \x\ < a 



'Zn c iz„a c iz n (x-a) 

Z n ~\~Z n 



(2.11) 



(2.12) 



, x > a 

for even n, where B was set to 1/2. This can be seen by using 

c i2~z n g Zn-Zn = f+1 , n odd 

z n + z n 1—1 ,neven, 

which directly follows from equation (2.8) with —inly; added to the exponent on the left hand 
side. 

We conclude this chapter with Figure 2.2, which illustrates that the Gamow functions do indeed 
increase exponentially in x. In the next chapter we will establish a connection between these 
functions and the generalized eigenfunctions. For this purpose, the explicit expression for the 
resonances will become relevant again, since we will discover that the generalized eigenfunctions 
have poles located at the resonances. 



2.4 Resonances and Gamow Functions 
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Figure 2.2. Plot of the squared modulus |Gi| 2 of the Gamow function corre- 
sponding to the first resonance with parameter choice a = 2 and A = 30. Notice 
that for x > a the function is strongly compressed in the x-direction, to make 
the exponential increase better visible. This was done by inserting the factor A 
in front of (x — a). 



CHAPTER 3 



Generalized Eigenfunctions and Strategy 



In this chapter we will establish precisely how the Gamow functions are connected with the 
generalized eigenfunctions of the potential V\. Thereby, the relation 

<f>(k,x) ~r](k)G(x) 

that was announced in the introduction will be turned into a rigorous statement, from which the 
function r\ will be deduced. But the major part of this chapter is devoted to the properties of r\ 
that are needed to prove the exponential decay of truncated Gamow functions within a certain 
time regime. 



1. Calculation of the Eigenfunctions 



First of all, the generalized eigenfunctions that belong to the potential V\ will be calculated. 
Clearly, there are symmetric and anti-symmetric solutions to a symmetrical potential and con- 
sidering both cases separately simplifies calculations considerably. An ansatz for the symmetric 
generalized eigenfunctions is given by 



<f>(k,x) 



Ae ikx + Be~ ikx 

Q{^ikx _|_ g— ikx^ 

Be ikx + AeT ikx 



, x < —a 
, \x\ < a 
, x > a, 



where k = y/k 2 — 2A. Due to the symmetry condition, it is sufficient to require that <f> and its 
first derivative are continuous at x equal to +a in order to guarantee continuity everywhere. And 
this yields the following homogeneous system of linear equations 



Ae 



-ika 



+ Be ika - C(e 



ika 



Ae -ika + B( ,ika _ k C ( e ika _ e ~i^ 

k 



-ika\ 




0. 



Upon solving this system, the symmetric generalized eigenfunctions are given by 

, x < —a 



4>(k, x) = Ae 



-ika 



p ik{x+a) _i_ F(k) —ik(x+a) 



F(k) COS kx 



„—ik(x—a) _i_ F(k) ik(x-a) 
e F{k) e 



Ixl < a 



. x > a 



(3.1) 



with ^(A:) = cos ka 



i— sin ka. 
k 
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Since this is a L°°-solution to the Schrodinger equation, it can not be normalized in the usual L 2 - 
sense. But similarly to the usual Fourier transform it is required 

that J 4>(k', x)4>(k, x) dx = 8{k! — k). And this gives A - 



4tt 

Analogously the ansatz for the anti-symmetric generalized eigenfunctions leads to the following 
expression 

ik(x+a) _ J(k) —ik(x+a) „, ^ __ 

J(fc) > X ^ " 



</>(k,x) 



4tt 



2i 



jjj^smkx ,\x\ < a (3.2) 

_ e _ ik(x _ a ) + m e ik(x-a) ; x > a 



k 

with J(k) = — cos ka — i sin ka. 

2. Refined Strategy 

The purpose of this section is to refine the argument that served to illustrate in the introduction 
how the exponential decay of e~ tHt tpo arises within a certain time regime. 

A vital part of this argument was the relation 4>(k, x) ~ r](k)G(x). But given that the generalized 
eigenfunctions solve (e — Hq)4> = V(f) their generic structure is 

<p{k,x) = e ikx + ${k,x), (3.3) 

that is general solution to the homogeneous part plus particular solution to the inhomogeneous 
part. And usually it is not the <£(/c, x) part that governs the time evolution 

e - iHt Mx) = J J r Mk)4>(k,x)e- i ^ t dk, (3.4) 

but the e lkx part. However, only $ can cause the generalized eigenfunctions to mimic the shape 
of the Gamow function G. Therefore, we are interested in the situation in which this function 
gives the dominating contribution to the generalized eigenfunctions. And heuristically this will 
be the case when 

Mk,x)\ = \4>{k,x)-e ikx \ > \e ikx \. 

Clearly, the generalized eigenfunctions that belong to V\ have the structure (3.3). And according 
to the above-mentioned argument, we expect that they take the shape of G when 

and |J| _1 »1 

respectively. In fact, the squared modulus of F _1 consists of a series of Breit-Wigner distributions 
whose height is of 0(A). This is illustrated in Figure 3.1 and will be proven rigorously later. Since 
these distributions are represented by the formula 

|«-i| 2 l a -i| 2 
(k — z' n ) 2 + z" 2 (k — z' n — iz'n){k — z' n + iz") 

their appearance in |-F| -2 reflects the fact that the continuation of (f>(k, x) to the complex fc-plane 
has poles due to F -1 . Therefore, we just need to apply the calculus of residues to the above- 
mentioned integral (3.4) to extract a term being proportional to exp(— z^z" t). So for those t 
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O(A) 




|F| 



-2 



1 



3tt 
2a 



57T 

2a 



Jfe 



7tt 
2a 



Figure 3.1. Plot of |F|~ 2 with parameter choice a = 2 and A = 100, where kF(k) 
denotes the Jost function of the symmetric generalized eigenfunctions. The fc-axis 
was rescaled by using the relation k = \fk 2 — 2A. 

for which this term gives the dominating contribution to e~ iHxt ij)o, the wave function will decay 
exponentially. But before we pursue this idea further, the existence of the poles will be proven 
in the next lemma. 

Lemma 4. The complex continuations of F and J have roots, which coincide with the reso- 
nances z n corresponding to odd n and even n respectively. In both cases the multiplicity of the 
roots is one. 

Proof. As already noted in Section 4, the resonances z n satisfy the following identity 



holds for odd n. And it is not difficult to see that this equation is equivalent to tan z n a = — zjp, 
which is the equation satisfied by the roots of F. An analogous argumentation holds for J with 
even n. 

To show that the multiplicity of the roots is one, we need to show that F'(z2 n -i) an d J'(z2n) do 
not vanish. In this regard, notice that 




, n odd 



n even. 



Using this, a straightforward calculation shows that 




F'(z2n-i) = —io, I cos ka — i - sin ka 



) 



, 2A 

l kk 2 



sin ka 



k=z 2n -i 
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Using the fact that F{z2 n -i) = 0, this can be expressed as 

n// \ / i\2\~ 
F {z 2n -i) = ~ [a - - ~- an ka 

\ fc / kk k=z 2n -i 

The term in brackets can not vanish, because Re Zm-\ is non-zero. Since the sine has roots only 
on the real axis and Im Zm-\ > 0, it can not vanish either. Hence, F'{z% n -\) is non-zero. An 
analogous argumentation holds for J. □ 

According to the previous lemma, the poles of the generalized eigenfunctions coincide with the 
resonances, which justifies that the Gamow 'eigenvalues' have been termed this way. In order 
to see this notice that kF(k) and kJ(k) are so called Jost functions (see [9, Section 2]), which 
define the Scattering-Matrix by 

q/ t \ F(k) J{k) 

S{k) = Fjk) and S{k)= J(k) 
respectively. Since it is common to define resonances as poles of the complex continuation of 
this operator, we see that they coincide with the roots of F and J and thereby with the Gamow 
'eigenvalues'. Although, this was only verified for V\, Skibsted's article [9] shows that it remains 
true for more general potentials. 

Another consequence of this observation, is that the Gamow function G n and the residue of 4>(k, •) 
at z n are linearly dependent. In order to see this let <f> be defined by 

= y(fc,Q 
> H(k) ' 

where H denotes either F or J. Then Res 4> = a_i tp(z n , •). And given that the Jost function 
evaluated at z n can be expressed as (see [9, Section 2]) 

G n (x) (p(z n ,x) 
£G n (x) i;<p(z n ,x) 

the linear dependence becomes immediately evident, when we recall that H{z n ) = 0. There- 
fore, we can establish a precise relation between the generalized eigenfunctions and the Gamow 
functions by expanding <p in a Laurent series about the resonance z n . This yields the expression 

4>{k, x) = — c G n (x) + Remainder, 

J£ Zn ^ 
V (k) 

from which r]{k) can be simply read off. 



3. Laurent Expansion of the Eigenfunctions 

This section will show that for big enough A there is an interval [z' n — 5,z' n + 5] on which the prin- 
cipal part r](k)G of the Laurent expansion of <j)(k, •) is a good approximation for the generalized 
eigenfunctions. In the next chapter this will allow us to justify the approximation 

e - iH ^ip Q = y^o(fc)0(^-)e~ i ^ t dA;~G(-) J ^VoC&M&K^' dk 

for those initial wave functions, whose energy distribution J-'tpo is concentrated in the above- 
mentioned interval. 



3.3 Laurent Expansion of the Eigenfunctions 



19 



Due to the fact that V\ has compact support, the generalized eigenfunctions break up in a part 
being valid on the support of V\ and a part being valid on its complement. Therefore, it is 
convenient to consider the time evolution integral (3.4) on these subsets of R separately. For 
now, we focus on supp V\. And this will require frequent use of the characteristic function x { _ a a] , 
which is why it is advantageous to introduce the following shorthand 

Moreover, it is convenient to restrict ourselves to symmetric initial wave functions such that only 
the symmetric generalized eigenfunctions (3.1) contribute to the integral (3.4). However, every 
result has a straightforward extension to the anti-symmetric case. 

Lemma 5. Let (f) denote the symmetric generalized eigenfunctions of H\ with X € M + and 
let Z2n-i with n > 1 denote its resonances. Then 4> can be expanded in the Laurent series 

<f>(k, x) = e~ ika ( u a - l{ - X) + ao(x) + R (k, x)) , (3.5) 

\K — Zm-\ J 

for every x in [—a, a]. And if X is chosen big enough, there is a 5 > 0, such that this expansion 
converges absolutely and uniformly on [z' 2n _i — 8,z' 2n _i + 5], provided n satisfies 



1 < n < J- - V2X =: n\. 



Proof. In Lemma 4 it was proven that the resonances Z2 n -i are first order poles of the 
complex continuation of <p(k, x) to the lower complex k plane. Thus, it is clear that the Laurent 
expansion of <fi is given by (3.5). In order to prove that it converges in the above-mentioned 
interval about z' 2n _i, we will determine its annulus of convergence. Since the resonances are 
separated points in the complex plane, it is clear that it will be a punctured disc around that 
particular resonance Z2n-i- So it remains to estimate the outer radius of convergence, which can 
be done by estimating the distance between neighboring poles of 0. 



By Lemma 4, the symmetric generalized eigenfunctions have poles numbered by odd indices, 
while the anti-symmetric ones have poles numbered by even indices. Therefore, we need to 
ensure that 

\Z2n+l ~ Z 2n ~\\ 2 > ±41+1 > ^4n-l f ° r ™ > 1, (3.6) 

where the last inequality will automatically hold if the first one is true. This is due to z^ 2 being 
a monotonically increasing function of k. Notice, that the factor in front of the imaginary part is 
chosen for convenience, but is not optimal. Using the explicit formula for the resonances given 
in Lemma 3, we find 



|^2n+l — Z2n-l\ 



(2n + l) 2 vr 2 



'2A + 



4a 2 



(2n- l) 2 vr 2 



2A + 



4a 2 



+ 



7T 



4a 2 2A 



2n + 1 



2n-l 



\ \ L + (2n+l) 2 7r2 ) (2n-l)27r2 ) 



1/2 



\ 2 
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By exploiting the formula (a + 6) (a — b) = a 2 — b 2 , the term in brackets appearing in the first 
summand can be simplified to 



TT 



8n 



4a 2 



> 



TT 



8n 



Aa 2 



The second summand can be simplified by using a similar argumentation, which yields the 
following lower bound 

/ \ 2 

' ' 1 4(2n + f) 2 



ILL 
4^? 2A 



1 



4a 2 2A 



1/2 



/ 



2 \ 2 i 
7T \ 1 



4a 2 



2A ox , (2n±l)^ ' 
/A + 4a2 



(2n+l) 2 7r 2 

Together these lower bounds imply that inequality (3.6) is full filled, if 



TT 

4a 2 



(An) 2 + ^-4(2n + if 



2A + 



(2n+l) 2 7r 2 
4a 2 



> Az 



112 

2n+l 



/vr 2 \ 2 1 4(2n+i; 

2A 2A + M0 



4 



Hence, inequality (3.6) holds if 



(4n) 2 > (2n+l) 2 ((2n+l) 2 -l). 
A 

Some simple rearrangements and estimates yield the following inequality 



> n 2 



1 n + 



1 

4' 



where equality would hold for 



2A 



1 ± 



V 



2A 



Clearly, n_ will approach zero in the limit of A tending to oo. Therefore, only n + is of relevance 
when A > 1 and n > 1. And if we regard n + as arithmetic mean, this allows us to conclude 
that the distance |^2n+i — Z2n-i| between neighboring resonances is larger than 2|z2 n+1 | for all 
integer n, which satisfy 



1 < n < 



So for these resonances there is a non-zero positive 5, such that the Laurent expansion (3.5) 
converges absolutely and uniformly for all k in [z' n — 5, z' n + 5]. □ 



That the region in which the Laurent expansion (3.5) converges ceases to intersect the real axis is 
not surprising. It reflects the fact that the influence of the resonances fades the larger their real 
part grows. From a physical point of view this can be anticipated, since initial wave functions 
with high energy (ipo,H\ipo) impinge on the 'barrier' with increased rate. Thus, their lifetime 
on top of the plateau is reduced. 

But more importantly, Lemma 5 shows that 

-ika a -i - I 1 



r] n (k) = e %ka - with a_i = Res 



k — Z2n-1 z 2n-l \y/7TF(k) 
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which verifies the Breit-Wigner shape of |r/ n | 2 as announced in the introduction. 

The estimates on the remainder of the Laurent expansion, will involve contour integrals around 
circles in the complex fc-plane. But to handle them some properties of the corresponding contour 
performed by k = y/k 2 — 2A are needed. And the relevant details will be established in the 
following lemma by considering the function y/c + de lip , where < d < c are real constants. 
Basically, this function will perform a contour which looks like an egg whose symmetry axis 
coincides with the real axis. This becomes plausible by rewriting the argument of the complex 
square root in polar form r{ip) exp(iS((p)). Then the effect of the square root on the exponential 
is to compress the circle in the direction of the imaginary axis, whereas its effect on the radius 
r(ip) shifts its center and will cause the circle to be compressed along the real axis. Since the 
strength of these effects depends on <p, the resulting contour will roughly look like an egg. For 
the purpose of estimating the remainder, the only information we need is in which region of 
the complex /c-plane the contour performed by k can be found. And the preceding illustration 
suggests that it is an annulus. The following lemma makes this precise. 

Lemma 6. Let < d < c be real and let ip lie in [0,2-7r). Then the complex contour yfc + de l( f, 
remains within the annulus S enclosed by the circles 

z ut = y/c + rout e l(p with r out = y/c - y/c- d and 

Zin = a/c + r in C %V With r in = Vc + d - y/c. 



PROOF. The assertion follows, if we show that the continuous function f(<p) defined 
by | y/c + d e 1 ? - y/c\ remains within the interval [rj n , r ou t[. This can be done by showing that / 
takes its only maximum at <p being equal to tt and its only minimum at ip being equal to zero. In 
order to do that, it is enough to show that f 2 is monotonically increasing on the interval (0, 7r), 
while it is monotonically decreasing on (ir,2ir). 

The monotonicity can be deduced from the sign of the first derivative of f 2 . But this will be 
easier if f 2 is rewritten first. By using the equality (a + b)(a — b) = a 2 — b 2 , we see that 

d 2 



\Vc + de il P + ^/c| 2 ' 
which has the following derivative with respect to <p 

d 2 



That means, f 2 is increasing whenever f 2 is decreasing and vice versa. In order to determine 
which of these cases occurs for which ip, we need to carry out the differentiation of f 2 symbolically 

d v f 2 = 2 (Re y/c + de^ + y/c) d^Re y/c + de^ 
+ 21m y/c + de^ d^lm y/c + de il ?. 

With the help of the formula \Jx + iy = 2" 1 / 2 ({yj x 2 + y 2 + x) 1 ^ 2 + isgn(y)(A/x 2 + y 2 — x) 1 / 2 ^ , 
which is easily proved by rewriting the argument of the square root in polar form, we can bring 
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real and imaginary part in the following form 

Re \/ c + d e lL f = — ^ \J c 2 + d 2 + 2cd cos ~tp + c + d cos </?^ ^ and 

Im v c + <ie^ = sgn(sm<^) — c 2 + <i 2 + 2cdcos — (c+ <icos 99) J 
Using these formulas, the real and imaginary part can be differentiated with respect to ip, yielding 

o r, / , — ~ 11/ —cdsmp , . \ , 

a w Re v c + o e* 1 ^ = . — . — d sin ip and 

4 Re y/c + de^ \ ^c 2 + d 2 + 2cdcosp J 

. T / , — — 1 sgn(sinco) / —cds'mp , . \ 

4 im y/c + de** \ yjc 2 + d 2 + 2cd cos <p I 



Therefore, we arrive at the following expression for the derivative of / 



2 



~ 2 ds'mp I c , 1 I / 1 , 



a,/' = _= _ - ^ = + 11 + 



y^ 2 + d 2 + 2cdcosv? /V Re y/c + de** 

ds'mip I —c \ . t . 

+ "^r~ 7 9 J9 ^ + 1 sgn(sm^) 
2 y ^c 2 + d 2 + 2cd cos 99 J 

= g(<p) + /i(y>). 

Comparing both summands factor by factor, we see that \g\ is strictly bigger than \h\ on the 
interval (0, it) and (71", 2ir). And the signs of g and /i are such that f 2 is monotonically increasing in 
the first case, while it is monotonically decreasing on (it, 2tt). Therefore, f 2 reaches its maximum 
at 7r and its minimum at zero. So if we define 



r ou t = /(tt) = y/c - \Jc— d and 

Tin = /(0) = Vc + d - y/c, 

we can conclude that the contour performed by Vc + d e tL P remains within the annulus S, which 
is enclosed by the circles with radii rt n and r ou t both being centered at \fc. □ 

We can now estimate the modulus of the remainder of the Laurent expansion. 

Theorem 1. Let 4> denote the symmetric generalized eigenfunctions of H\ and let n be a fixed odd 
integer satisfying 1 < n < n\. If X is chosen big enough, the Laurent expansion of e lka cf)(k, •) in k 
about the resonance z n converges for all k in [z' n — 5,z' n + 5] with 

<5= \{z n -z n _ x ). 

Furthermore, the Laurent coefficients a m (-) and the remainder Ro(k, •) satisfy 

||x«am(0lloo <0{\ m ' 2 ) and 
\\XaRo{k, ■) ||oo < 0{\ 1 ^ 2 )\k - z n \. 

Proof. The L°° bounds on the coefficients will be proven first. Since (3.5) is a Laurent 
expansion its coefficients can be calculated by 

If 6(z.-)e iza 
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where C is a closed contour around the pole z n , which remains within the annulus of convergence. 
Using the explicit form of the generalized eigenfunctions given in (3.1), we therefore have the 
following upper bound 

11 

||XaOm(-)lloo < —F= SUp z T ~—. — SUp SUp \Xa COS Zx\ T m . (3.7) 

V 71 " zee cos az - i- sin az\ z &c xeM 



Here C was chosen to be the circle 

{z n + re iv \<p <E [0,2?r)} with (3.8) 



z. 



j _ i 



which is a contour within the annulus of convergence for all n, since only one pole is enclosed. 
To handle the suprema appearing in (3.7), we need to know how the circle C in the z plane looks 
like in the z plane. Due to the fact that we are only interested in the limit of A tending to oo, 
we can exploit the dependence of z n on A in order to find 



z 



(z 2 - 2A) 1/2 = (4 - 2A + 2z n re i{p + r V 2 ^) 1/2 
= (# + 2« + i<)re*) 1 / a 



, 2 / 1 2n-l 2 vr 2 1 2n-2 2 vr 2 \ , ^ 1/2 

z 2 + 2z' ri V2X 1 + ~i =-4 1 - > ^ 

2 4a 2 2A 2 4a 2 2A 



which does not depend on A anymore. Thus, 



sup sup \xa cos zx\ < sup sup Xa cosh z" x = sup cosh z"a < c, (3-9) 

z^Cxm z&c iei zee 

since z" is independent of A on the circle C. And for similar reasons 

H|sinaz| < ^ =0(A- 1 / 2 ), 
\z\ \z\ 

from which we can deduce 

l (. ,~z . A" 1 / .. 1 



sup = < inf | cos az\ — | - sin az\ \ < inf | cos az l 

z&c | cos az — i* sin az\ \zec z J \z6C 4 

provided A is big enough. So it remains to estimate the infimum of |cosaz| on the circle C. 
Regarding this, the information given in Lemma 6 will be sufficient, due to the maximum modulus 
principle of complex analysis (see e.g. [6]). According to this principle, the modulus of an analytic 
non-constant function considered on a connected open subset S of C, reaches its maximum on 
the boundary dS. Lemma 6 proves that the contour C remains within the annulus S bounded 
by the circles 



Zin = z n + ( \/2z'^ - z®_i - z n ) e lip and 



Zout — Z n + {z n Z n -\) e ^ • 

Since these circles enclose exactly one root of cos az at their center, (cos az)" 1 is analytic within 
the annulus S. Thus the maximum modulus principle is applicable and it can be concluded 
that I cos azl reaches its minimum on one of the circles Z{ n or z 0U f. 
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In order to find the minimum of | cos az\ on the above mentioned circles, it is easiest to consider 
the squared modulus. Since az' n is a root of the cosine, the squared modulus can be rewritten as 
follows 

| cos(az' n + are lip )\ 2 = \ sm(are tLp )\ 2 = - (cosh(2ar sin if) — cos(2ar cos if)), 

where r denotes either rj n or r out . Due to the symmetry of this expression about w, it is enough 
to consider 

dip\ s\n{are llp )\ 2 = ar(sinh(2ar sin if) cos if — sin(2arcos if) sin if) 

on the interval [0,7r]. For if in (0,7r/2) we also know that cos if and thereby the first summand 
is positive. In this case 

d v \ sin(are i</3 )| 2 > ar sin if(2ar cos if — sin(2ar cos if)) > 0, 

since sinhx is bigger than x and (x — sinx) is bigger than zero for all positive x. However, 
should ip lie in (tt/2,tt) we know that cos if is negative. Then 

dp\ s\n{are llp )\ 2 = — ar(sinh(— 2ars'mf) cos if — sin(— 2ar cos if) sinif) 
< — ar sin if((—2ar cos if) — sin(— 2ar cos ip)) < 

for the same reasons as above. The monotonicity implied by these considerations shows 
that | sin(are J</3 )| 2 reaches its minima at zero and it. But due to symmetry the modulus of this 
function at zero will not differ from its modulus at tt, which is why 

inf |cosaz| > inf |cosaz| = min{| coso5j n (0)|, | cos az ou t(0)\} ■ 

z(zC z£S 

A lower bound for this minimum can be found by determining lower bounds for the cosines. 
Since Zi n (0) = ^2z'^ — z l2 _ x , we see that 



cos(az m (0))| 2 = ^1 + cos (2a^2~z% -~z' 2 _^j = * (l + cos Umr^Jl^ 



2n 2 



where the explicit location of the resonances z n was used. Considering the square root in the ar- 
gument of the remaining cosine as deviation from 2nir, where the cosine would take its maximum, 
it becomes clear that 

| cos(az m (0))| 2 n > *(1 + cos(V2vr)) > l - (l + cos (^j ) = \ (\ - 



Using the explicit location of the resonances z n again, we also find 

| cos(az out (0))| 2 = | cos(n7r)| 2 = 1. 
And this immediately implies that 



inf I cos az\ > 




from which we can conclude 



sup ^r—^ — < ( inf | cos az\ — - ) < 12 



zee | cos az — i| sin az\ \zec I 
provided A is big enough. 
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This was the missing estimate for the upper bound of the L°°-norm of the Laurent coeffi- 
cients XaO"m(')- Using this inequality together with inequality (3.9), we can deduce from the 
usual Cauchy estimate (3.7) for the Laurent coefficients that 

IIXaOmOHoo < CT~ m = C (z' n - Z^)'™ = 0(X^ 2 ). 



It remains to prove the bound on the remainder Rq of the Laurent expansion. Multiplying (3.5) 
by e lka and subtracting the principal part of the Laurent expansion on both sides yields 



e ika ^{k,x) - 



a-i(x) 



= ao(x) + Ro(k, x). 



k z n 

Since the right hand side is a Taylor expansion of the left hand side in the complex variable k, 
the Lagrange form of the remainder can be used, to get an upper bound on the modulus of Rq. 
(Notice that, the Lagrange form of the remainder only gives a bound in dimensions bigger than 
one, but it is not exact anymore.) So if k is within the radius of convergence of the Taylor series, 
there is a 



C £ {z n + t(k — z n )\t 6 (0, 1)} such that 



Xaflo(MHoo< Xad k [e tka <P(k, ■)](() 



\k- 



where we exploited that the a m are equal to the Laurent coefficients of e lka (j)(k, •). Using Cauchy's 
formula, the derivative that appears in the above bound can be estimated as follows 



X<A[e'*V(*,-)](C)Hc 



1 

2n 



<j iza 4>{z, •) 



dz 



< — supsup\x a e lza <j)(z,x)\ 
^ zee xm 



1 



dz, 



lc\z-C\ 2 

where C needs to be a closed contour that contains Q. The supremum can be estimated as above, 
when the domain of £ is contained in the interior of the circle (3.8), provided C is chosen to be 
this contour. If we postpone proving that Q satisfies this condition, we can conclude that the 
supremum of \\xa,e lza '<fi(z , -)||oo over all z £ C is of 0(A°). And this together with the particular 
form of the contour C, implies 

||x<A[e*V(£v)](C)Hoc < — ^ d<p 



2ir 



< O(A 



inf 

V6[0,2tt) 



\z n -C + re^\ 2 



where r denotes (z' n — z' n _^). The term in brackets can be estimated as follows 
\z n ~ ( + re i<p \ 2 = \w + re iip \ 2 = (w' + r cos cp) 2 + (w" + r sin tp) 2 

> (| w '|-r) 2 + (|«/'|-»") 2 
>2r(r-(\w'\ + \w"\)). 

Since k is in [z' n — 8, z' n + S] with 5 being equal to \{z' n — ^_x), the real part of ( will lie in the 



very same interval, which implies 

\w + re^\ 2 >r{z' n -z' n _ l 
Thus, we can conclude that 

||Xoi2o(Mlloo <0{\ l/2 )\k 



2\w"\) > r(r-2\4\) =rO(\- 1/2 ). 
n \ Vk£[z' n -5,z' n + 5], 
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provided A is big enough. So it remains to show that the domain of ( is contained in the interior 
of the circle defined by (3.8). This is the case if 

\k - z n \ 2 < (z' n - 4„ x ) 2 Vk g [z' n - 5, z' n + 6}. 

But on this interval 

\k - z n \ <o + z n = -r + z n . 
So for fixed n, the variable C will lie in the interior of the circle (3.8), if 

< 2 <^ 2 = 0(A- 1 ). 

For big enough A this condition is full filled, since z" 2 is of 0(A -2 ). Notice that r is smaller 
than \z n — z n _2\- Hence, this condition will also ensure that the k interval of interest, will lie 
within the annulus on which the Laurent expansion of 4> converges. □ 



CHAPTER 4 



Exponential Decay 



The purpose of this chapter is to show that the truncated Gamow functions XaG n decay expo- 
nentially within a certain time regime on the support of V\ . Thereby, the last step of the heuristic 
discussion given in the introduction, which is 

e- iHt ip (x) ~ J J r Mk)v(k)G(x)e- i ^ t dk ~ G{x) e _ H (4.1) 

will be turned into a precise statement. We will discover that the error which is produced by this 
approximation decreases with increasing A. And after we have determined the regime in which 
the exponential law is valid, we will conclude by extending these results to initial wave functions 
which lie in C 2 ([—a,a\) and vanish at ±a. 

1. Truncated Gamow Functions 

The exponential decay of the truncated Gamow functions XaG n within a certain time regime 
follows from the exponential decay of any other initial wave function lying in their neighbor- 
hood (with respect to the L 2 -norm). This is due to the unitarity of e~ iHxt , which implies that 
an L 2 -neighbor ip n of the truncated Gamow function XaG n satisfies 

We'^XaGn - e-^Vn||2 = || XaG n - V>n|| 2 < 1. 

So for now it will make no difference if we choose the bound states of H\ restricted 
to Ld = {Y> € L 2 {[— a, a])\ip(a) = 4>{—a) = 0} in place of XaG n . But we will benefit from 
this choice when the results derived for these specific initial wave functions will be extended to 
generic ones. Apart from that, we restrict ourselves to symmetric bound states for convenience, 
and in order to see that they lie in the vicinity of XaG n , notice that the symmetric bound states 
are given by 

Y> n = C0S(4„-1 -)Xo- 

Therefore, they differ from the symmetric truncated Gamow functions only by the appearance 
of the real part z' 2n _i instead of the entire resonance Zm-\- Given that the imaginary part of 
the resonance is of 0{\~ 1 / 2 ) we only need to choose A big enough to ensure that the distance 
between the symmetric bound states and the symmetric truncated Gamow functions is small 
enough with respect to the L 2 -norm. 

In order to show that the time evolved symmetric bound states decay exponentially within a 
certain time regime, we need to determine e~ lHxt tp n - In this regard, we will exploit that for 
big A the major part of the L 2 -mass of the energy distribution Fip^k) is concentrated in a small 
interval about z' 2n _ 1 . This becomes evident from 

, / , \ o , , v Cn G- COS JvQj , , „ , 

F ( k ) Z 2n-l- k2 
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Its second factor resembles a sinc-function which is peaked about z' 2n _ x an< ^ nas ro °ts at z' 2m _i 
for all m^n. Due to the fact that \F(z' 2m )\- 1 is of O(A ) and |F(4 m _ 1 )r 1 is of C^A 1 / 2 ) for all 
integer m, this causes |V>„| to be strongly peaked in a neighborhood of z' 2n _i while being of 0(A°) 
everywhere else. And since the distance between neighboring roots of the energy distribution is 
of 0(A -1 / 2 ), we conclude that the major part of the L 2 -mass of V> n is concentrated in the interval 

[4n-i-**2„-i + <*] with 5= l -{z' 2n „ l -z' 2n _ 2 ) = 0{\- 1 ! 2 ). (4.3) 

This property of the energy distribution ip n is the reason why the time evolution of the symmetric 
bound states is governed by the generalized eigenfunctions 4>(k, •) that belong to the interval (4.3). 
So this property is the reason for (4.1) to be a good approximation. However, for the rigorous 
justification of (4.1) we need to estimate the error that is produced when (f>(k,-) is replaced 
by Vn(k)G n , which is its leading order contribution on the interval (4.3). Therefore, we will 
consider 

r°° k 2 

Xae~ iH ^ n = Xa / M^Vni^Gn^e'^ dk 
Jo 

+ Xa / Mk)Wk, ■) - r/n(fc)G n (-))e"^' dk 
Jo 

and prove that the L 2 -norm of the error integral 

r°° - fe 2 

E n (-,t) = X a / Mk)(<P(k,-) ~ Vni^Gni-^e-^ dk 
Jo 

is small, while we will see that its main contribution 

*n(", t) = Xa / Vn(fc)r?n(fe)G„(-)e-^ t dk 

Jo 

decays exponentially within a certain time regime. 



2. Truncated Gamow Functions - Main Contribution 

In the current section it will be shown that the main contribution ^ n to Xa^ lHxtr ^n decays 
exponentially within a certain time regime and moreover this regime will be determined. 

In order to prove that there is a time interval in which 

*n(-, t) = XaG n / Mk) e~ ika 1 e"^ dk 

Jo k - Z2n-\ 

decays exponentially, we will use the calculus of residues. Due to the term exp(— i^-t), the 
integration contour needs to run through the lower complex /c-plane in order to get converging 
integrals. In particular, we will choose the contour illustrated in Figure 4.1 with R tending 
to oo, which results in an exponentially decaying contribution to caused by the pole of r\ n . 
Therefore, our main task will be to determine the contribution from the rest of the integration 
contour and to show that it can be neglected during a certain time regime. 

Lemma 7. Let n be a fixed odd integer satisfying 1 < n < n\. Then for all t > 0, we have 

*n(-, t) = cXoG„e-4* + Xa G n r(t), (4.4) 



4.2 Truncated Gamow Functions - Main Contribution 
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Figure 4.1. Illustration of the integration contour C in the complex fc-plane, 
where R is understood to tend to oo. 

where the non-zero constant c is of O(X ) and 

\r(t)\ < Oix-^r 1 / 2 . 



Proof. By the residue theorem 



c 



Tp n {z)r] n {z)e dz = -2m Res (i/i n r/ n e l l^ 2 ') 



where C is the contour illustrated in Figure 4.1. From this figure it becomes clear that the 
contour integral breaks up into 



f 

Jo 



2 

TpnityVnitye'^* dk + 



C1+C2 



. z 2 

*l>n(z)Vn(z)e~ %li ~ t dz. 



We will show first that the arc-like part C\ of the integration contour gives no contribution. In 
this respect, notice that 

I / I < lim R ^{Re-^niRe-^e-^ 2 sin2 <^. 

J Cl R^oo J 

As lim^oo \i\) n {Re~ %ip )e~ %aRe lv \ = for all ip in [0,7r/4], which becomes evident by a straight- 
forward calculation, we only need to choose R big enough to ensure the following estimates 

I f I < lim cR r / \(Rcos<p- z'J 2 + (Rsinip + Z y)- 1 / 2 e- t 2 R2sin2 ^dip 

J Cl R^oo J 

/■7T/4 t 22 

< lim cR / |i?cos (p - — cos ^j" 1 e~i R ^ dip 
R-^oo J 2 



< lim c' / e~2 H * <p d<p = hm c" 



R^oo 



R^oo 



tR 2 



0. 



Let us estimate the modulus of the contribution from C2 next. Using the explicit form of ip n 
and r] n , the following estimates are immediately evident 

\c n \ 1 



c 2 



< 



< 



\c n a. 



yl2 



H 



|1 + i- z tanaz| \z'n + ir z \ \z — z n \ 

00 



a_i _t 2 

e 2 dr 



1 + z-tanaz z — z n \ e 2 dr. 
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To handle the first factor of the remaining integrand, we will show that Re (i| tan az) is positive, 
since this implies the following inequality 



i-i - z \ z ~i-i 

\1 + i- t&n az\ < -tanaz 

1 z z 



In order to see this, notice that 



a rg (if tan d) = + a rg ( !V ^T^) + a rg (tan al) and 



sin2az +isinh2az 

tanaz = — . 

cos 2az' + cosh 2az" 

Therefore, z" > y/2X and this implies that 

arg(tanaz) G ^7r^ as well as | tanazj > 0(A°) / 0. 

But then arg(z| tan az) remains within the interval (Stt/2, 2ir + 7r/4), which justifies 

/ < c « — / \ re z n\ e 2 r 

■/c a 4 Jo 

provided A is big enough. By determining the maximum of the first term in the remaining 
integrand and using Theorem 1, we can thereby conclude that 



JC 2 Z n \ z n ~ \ z n\) JO 



So it remains to determine the residue and its order with respect to A. First of all, it is not 
difficult to see that 



Res(ip n r) n e ) = a_i^ n ( 

Zn \ / 



2 2 



From the previous section it is known that IVVi^n)! ^ s °f ^(A 1 / 2 ). Since z^ is of 0(A _1 ), the 
continuity of ijj n in the region between z n and the real axis requires |?/> n (z n )| to be of 0(X 1 ^ 2 ). 
But according to Theorem 1, we have |a_i| < 0(\~ 1 / 2 ). So we conclude that c is of 0(A°). □ 

To determine the regime in which \l/ n decays exponentially consider the ratio of the two summands 
that contribute to \l/ n (equation (4.4)), that is 

|r(f)| 



z 2 

— 4 

ce 2 



< 0(A" 5 / 4 ) e r "* (Tnt)- 1 ' 2 . 



According to Lemma 7 exponential decay prevails, if this quotient is much smaller than one. 
Clearly, this can not be the case for t <C 1 nor for f > 1. But this is exactly what we expected, 
since it was shown in Section 2 that the survival probability |(Y>o, e ~ tHt i ; o}\ 2 can n °t decay 
exponentially in these time regimes as long as tpQ is a non-zero element of D(H). And the initial 
wave functions ip n are elements of V(H\), because their weak derivatives tp' n and exist and 
have finite L 2 -norm. However, on intermediate time scales the above-mentioned quotient is small 
for periods of several lifetimes provided A>1. This becomes immediately evident when t is such 
that it satisfies 

1 < T n t < N 
with an integer N > 1 that depends on A. 



4.3 Truncated Gamow Functions - Error 
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3. Truncated Gamow Functions - Error 



To finish the proof that e lHxt ip n decays exponentially in a certain time interval on the support 
of V\, it remains to estimate the L 2 -norm of the error integral 

E n (;t) = X a / Mk)Wk,-) ~ Vn^Gr^e-^ dk. 

Jo 

We will see that this norm is of 0(A -1 / 4 ). Therefore, the main contribution determined in the 
previous section will approximate Xa&~ lHxtr ipn better the bigger A is chosen. This is in accord 
with our intuition that an initially localized wave function remains longer within the support of 
V\ the higher the potential plateau is. 

Let us illustrate briefly why we should expect E n (-,t) to have small L 2 -norm before it will be 
proven rigorously. It was explained in Section 4.1 that the energy distribution tp n is strongly 
localized within the interval I given by [z^n-i ~~ ^ z 2n-i + By exploiting this fact we will 
be able to show that the contribution from M + \ I to E n has small L 2 -norm. To estimate the 
contribution from X itself, we will use the fact that r] n G n is the principal part of the Laurent 
expansion of (j). But since this part of the proof is a bit more involved, we will briefly outline the 
approach. Using Lemma 5 and Theorem 1, the contribution from X can be expressed as 

Xaao(-) / Mtye-^e-^dk + xa [ Mk) e - lka R (k, .Je"*^ dk. 
Ji Ji 

The details on the remainder Rq proven in Theorem 1 are sufficient to handle the second integral 
in a straightforward manner. The involved part is the first integral. In this respect, notice that ip n 
contains F -1 whose meromorphic continuation has poles at z n instead of z n . Therefore, F _1 has 
a Laurent expansion which is completely analogous to the one of F _1 , namely 

F-\k) = -^ + R(k). 

K Z n 

From this expression it becomes clear that the phase of the principal part of ip n oscillates rapidly 
on the interval I. Therefore, we will use a stationary phase argument to handle the major 
contribution to the first integral from above. 

Lemma 8. Let n be a fixed odd integer satisfying 1 < n < n\. Then for all t > 0, we have 

\\E n (;t)\\2<0(X-^ 4 ). 



PROOF. For brevity let I denote [z' n — 5, z' n + 5], where 5 can be chosen as \{z' n — z' n _^) due 
to Theorem 1. Then, 



\E n (;t)h < 



Xa (1 - Xx )(k)Mk)m, •) - TM(*0G„(-))e-*Tr* dk 
Jo 



Xa ^ Mk)(Hk, ■) ~ %(fc)G n (-))e-^* dk 
A + B. 
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First consider A. The following estimate is immediately evident 



A < 



poo 

1 1 



/•oo 

+ \\XaG n \\ 2 / {l- Xl ){k)\Mk)Vn{k)\dk 
Jo 

=■■ A x + A 2 , 

where Xa was dropped in the first summand. From the spectral theorem, we deduce 



-M x t 



(l- Xx )(k)Mk)<f>(k,-)dk 



-Mxt -t— 1 



Using the unitarity of T and e this reduces to 



(1 - Xx )(k)\Mk)r dk < / (1- Xl )(k) 



yl2 



dk < o(\- 1/2 ), 



k 2 y 



where the last step follows from a straightforward calculation. Regarding A 2 , notice that 

HXaGnlloo = ||XaCOs(5 n -)||oo < 0(X°), 

which follows from the A-dependence of z n given in Lemma 3. Therefore, we find the following 
estimates 



k£R+\I \ K ~ Z n\ JO \ z n — K I 



dk 



<c'/ (l-Xx)Wi 



dk 



|y'2 1-21 
10 \ Z n ~ K I 

< o{\- 1 / 2 ), 

where the last step follows again from a straightforward calculation. And this directly implies 
that 

A<0{\- l/A ) Vt>0. 
So it remains to estimate B. From Lemma 5 and Theorem 1, it is known that 



<f>{k, ■) - r, n (k)G n (-) 



-ika 



a (-) + Ro(k,-)) Vk£l 



with Hxaaolloo < O(A ) and \\ Xa Ro(k, OIU < 0(X^ 2 )\k - z n \. 
Thus, using the fact that x takes only values within a compact set we find 



B < O(A ) 
If we define 



4>n(k)e 



-ika e -i^t dk 



+ 0{\ 1 ' 2 ) 5 sup \4> n (k)(k - z n )\ =: B 1 + B 2 . 
hex 



h n (k) 



cos ka 



r /2 _ h2 

for brevity and apply the results of Theorem 1 to F _1 , we find 

sup \i> n (k)(k - z n )\ = \c n \ sup \F~ l (k)(k - z n )h n (k)\ 
k&x fceX 

< c' n sup I F-\k){k-z n ) | 
fee J 



< c' n sup(|a_i| + \a (k - z n )\ + \R (k)(k - z n )\) = 0{\- l/2 ) 
fcex 
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And due to the fact that S is of 0(A 1//2 ), we conclude 

As was described above, we want to apply a stationary phase argument to estimate B\. For this 
purpose we will use Theorem 1 once again in order to expand F^ 1 , which is part of tp n . This 
immediately implies 



B x < 0(A°) 
< 0(A°) 



h n (k)e 1 2 1 dk 



+ O(A ) 5 sup \{a + R {k))h n (k)\ 
feez 



fj n (k)h n (k)e 1 2 t dk 



+ 0(A 



-l/2\ 



where % was defined as a_i(/c — z n ). Since h n is real valued, it does not contribute to the phase 
of the remaining integrand. Therefore, we can use partial integration in order to get 



< 



\rj n (k)\h n (k)e iS W- i ^ t dk 



\f]n{k)\h n {k) lS (k)~i4 
d k S{k) - kt 



■ 4+<5 








- <S 





\fjn(k)\h n (k) 



d k S{k) - kt 



dk, 



(4.5) 



where S denotes the phase of fj n . Due to the fact that fj n is of 0(A°) at the boundaries of 
integration, as well as 



dkS(k) = dk arctan 



<0, 



h — yl (h — 7 M 2 4- y" 2 

«, z n z n ) -f z n 

we immediately see that the first summand in (4.5) is bounded by a constant of 0{X^ 1 ^ 2 ) for 
all t > 0. 

The second term that appears in (4.5) will be of the same order, if the supremum of its integrand 
is bounded by a constant of O(A ). An application of the product rule yields the following upper 
bound for the integrand 



X 



hn | Vn I 



T S'-(-)t 



+ 



X 



\Vn\K 



+ 



x 2 



S"-t 



I V 'n I ^T. 



Straightforward calculations justify the following estimates 

\a-i\(k - z' n ) 



Xj 



and ||x x <S"| 



X 2 



{{k-z' n f + <«)3/2 



< 0{\ 1 ' 2 ) 

> 

< 0{X 1 ' 2 ). 



{{k-z' n y + z'>^\ 

To get an upper bound for ||Xz^4lloo notice that h n is an element of the Sobolev space H 2 {E). 
We can therefore apply Sobolev's inequality (see [7, Chapter 8]) in order to find that 

RM*)II« ^ \ (H^M*)lll + \\d 2 k h n {k)\\ 2 ) < 1 (\\xMx)\\l + IkVn(x)Hl) < 0(1), 

where we used the fact that h n is the usual Fourier transform of ip n . Putting all these inequalities 
together we find 

I fjn | h n 



Xi 



S> - (-)t 



< 0(\ 1 ' 2 ) 

<o(i) 



S' - (•)« 
Vt > 0. 



+ 0(A) 



x 2 



(5' - 



+ 0(A b 



X 



x S*-{-)t 
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This shows that B\ has to be bounded by a constant of 0(X 1 / 2 ), which allows us to conclude 
that 

||£n(-,*)||2< A + BKOiX' 1 ^) Vt>0. 

□ 



4. Generic Initial Wave Functions 

From a physical point of view it is not satisfactory, that only the specific initial wave functions 
considered in the previous sections show exponential decay in certain time intervals. Therefore, 
we will extend the results achieved so far to generic ipo in this section. In this process we 
will understand why truncated Gamow functions are distinguished initial wave functions. And 
moreover, it will become clear why the bound states of H\ restricted to are better adapted 
for the extension to generic initial wave functions than the truncated Gamow functions. 

But what does generic mean in the first place? The probability that the meta-stable particle 
already decayed should initially still be zero. So we still demand that supp ipo C [—a, a] . For sim- 
plicity we will restrict ourselves to symmetric wave functions again. And since we exploited the 
strong localization of Ftpo, the Riemann-Lebesgue Lemma suggests that ipo should be sufficiently 
smooth. Therefore, under generic we shall understand that 



^0 G {/ G C^[-a,a])\f(x) = f(-x)} =: A k , 



where k will be chosen big enough and Cq([— a, a]) denotes the space of fc-times continuously 
differentiable functions on [—a, a] which vanish at ±a. 

The fact that the generalized eigenfunction <p(k, ■) could be replaced in (4.1) by the principal part 
of its Laurent expansion about z n without producing a huge error, was possible only since the 
contribution ||xiR+\[<-<5,<+5]-^ r ' ! / ; o||2 from all the other resonances was small. So for the method 
from the previous section to be successful, this one-resonance-approximation has to remain valid 
for elements of Ak- However, it is not difficult to see that this will not happen without any 
further restrictions on ipo, even when k = oo. Choose for example 

i/)q(x) = (cos z[x + cos z' 3 x)xa{x), 

then this will be an element of Ak and we know from Section 4.1 that 

_ , . e lka ( ci cos ka C3 cos ka \ 

But this is strongly localized about z[ and z' 3 , which implies that ||xiR+\[zj-(5,z^+<5]^ r V'o||2 is not 
small. This shows that the initial wave functions 

B = {cos(z' n -)Xa} nen 

considered in the previous sections are distinguished insofar as they allow for a one-resonance- 
approximation. And since the truncated Gamow function cos(z n ■)%„ belongs to the neighbor- 
hood of cos(z' n -)xa (measured in the L 2 -norm) it is distinguished for the same reason. But for 
generic initial wave functions there will be contributions from more than one resonance. And in 
order to handle them, we will follow the hint that is actually given by the above counterexample. 



4.4 Generic Initial Wave Functions 
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As already noted in Section 4.1, the initial wave functions Bq are a subset of 

B = i sm 5 f 

I 2a J neN 

which are the eigenfunctions of Hq with domain Cq([— a, a]). From the theory of Fourier se- 
ries (see [1, Example 1.2.3]) it is well-known that B is a complete orthonormal set in L 2 ([—a, a]). 
One consequence of this is the essential self-adjointness of Hq on its domain. But more impor- 
tantly we can conclude that Bq is an orthonormal basis in the set of general initial wave functions 
Ak- The counterexample given above is therefore just one particular element of Ak expressed as 
a Fourier series, which as we now know can be done with any other element as well. 

Thus, we can now determine the time evolution of our general class of initial wave functions 
through a term by term application of the one-resonance-approximation, which yields 

oo 

e~ lHxt ifj (x) = ^2a n C n (t)G n (x) + R(x,t) for all x G [-a, a], 

n=l 

where C n (t) = / J : ijj n {k)r] n {k)e'" l ^ t dk with ip n £ Bq. (4.6) 
Jo 

The exponential decay within a certain time regime then immediately follows from Lemma 7. 
And to avoid questions of summability we will exploit the smoothness of the elements in Ak ■ This 
property together with the Riemann-Lebesgue Lemma ensure that the modulus of the Fourier 
coefficients a n is as small as we please if only n is big enough. Therefore, it will be sufficient to 
apply the one-resonance-approximation to the first few summands, which represent the major 
contribution to i/jq. 

Before we will use the preceding considerations to prove the next Lemma, notice that contrary 
to Bq it is not evident that the set of truncated Gamow functions {xaG n } n&N is a basis in Ak- 
Therefore, it is not clear which set of initial wave functions actually is covered by finite linear 
combinations of them. 

Lemma 9. Let ipQ be an element of A\ and let C n {t) be defined as in equation (4.6). If X is big 
enough, there is an integer 1 < N < n\ such that for every x £ [—a, a] 

N 

e- iH ^o(x) = Y j a n C n (t)G n (x) + R(x,t) 
n=l 

with ||xo-R(-,t)||2 < 0(A~ 1/4 ) for allt> 0. 

Proof. Let ip n with n > 1 denote the elements of Bq. Following the considerations from 
above these elements build a complete orthonormal set in Ak- Hence, for every ipQ in Ak we have 

oo N oo 

e- lHxt ^Q = e- iHxt Y J {^^n)^n = e- iHxt Y J ^n + e-^ t £ a n ^ n 

n=l n=l n=JV+l 

= e- iH ^ + R(-,t), 
where the integer iV > 1 will be specified later. 

The Fourier coefficients a n are independent of A, since neither ipQ nor any of the ifj n depends 
on this parameter. Therefore, we can use the Riemann-Lebesgue Lemma together with the 
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differentiability of ipo, in order to see that for any m < k there is a A-independent constant C m > 
such that 

l i ^ C m w ^ 1 

\a n < V n > 1. 

This inequality together with the unitarity of e~ lHxt imply 



W;t)\\l 



Y a n^n = Y I ™' 2 " C k Y 11 
n=N+l 2 n=N+l n=N+l 

where Plancherel's identity was used in the second step. For this sum to converge, it is sufficient 
to choose k = 1. In this case we find the following upper bound 



N 



\\R{-Ml<cl Y / m ~ 2dm 

n=N+l Jn-l 

By Lemma 8 we can therefore find a 1 < N < n\ such that the one-resonance-approximation 
applies to every summand of e~ lHxt 4>o, while ||i?(-,t)||2 < 0(A _1 / 4 ) for all t > 0. 

And using Lemma 8 once again, we conclude that for these N and every x G [—a, a] 

N 



e- lHxt i; = Y a nMt)+R'(-,t), 



n=l 



where ip n (t) denotes the main contribution from the one-resonance-approximation applied to e ^^ x ^ipn 
that is 

.■k a 



roo 

Mt) = G n / FMVVnWe-^ dk 
Jo 

and \\R'(;t)\\ 2 < 0(A -1 / 4 ) for all t > 0. 



□ 



CHAPTER 5 



Outlook 



In the previous chapters we studied the physical significance of Gamow functions by analyzing 
the Schrodinger evolution of initially localized wave functions on the support of the potential V\. 
The piece missing to complete picture is therefore the evolution of these initial wave functions 
outside of the support. In this regard, Skibsted's results [9] suggest that \e~ lHxt i/jo(x)\ increases 
exponentially in x up to a radius R(t) which moves outward, since he proved 

e- iHxt Xa G n £ e- ltz '/\ R{t) G n with R(t) = a + Re(z n ) t. 

On the one hand this is an astonishing assertion, but on the other hand large radii R(t) correspond 
to early escape. Thus, the exponential decay in time on the support of the potential should cause 
an exponential increase in x, otherwise the continuity equation would be violated. However, 
Skibsted's result was proven in the L 2 -norm, which raises the question whether the pointwise 
exponential increase of the wave function is actually true. Proving this with the help of an 
expansion in generalized eigenfunctions seems to be more delicate than proving the exponential 
decay within a certain time regime on the support of V\. This is due to the fact that the main 
contribution to the integral 

e~ lHxt Mx) = J ^Mk)<P(k,x)e- i ^ t dk, 

which is its residue, is of 0(\~ 1 / 2 ) for x > a. Since the lifetime T" 1 increases with increasing A, 
this is sensible from a physical point of view. But from a mathematical viewpoint the problem 
becomes more demanding, because this is even a lower order than some of the error integrals 
obtained when the time evolution on the support of the potential was studied. 

The question whether the Gamow functions have physical relevance, which was the subject 
this thesis, is part of the larger endeavor to study the Time-Energy Uncertainty Relation. This 
question appeared due to the connection between energy width and lifetime of meta-stable states. 
And although, the a-decay is a well-established application of the Time-Energy Uncertainty 
Relation, this inequality is surrounded by many mathematical questions. One of them is its 
rigorous derivation, which is problematic given that no self-adjoint time operator exists. However, 
this as well as several other questions concerning the Time-Energy Uncertainty Relation remain 
as a challenge for future work. 
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